Distribution of the quantum mechanical time-delay matrix for a 

chaotic cavity 



P. W. Brouwer 

Lyman Laboratory of Physics, Harvard University, Cambridge, MA 02138, USA 



K. M. Frahm 

Laboratoire de Physique Quantique, UMR 5626 du CNRS, Universite Paul Sabatier, 31062 

Toulouse Cedex 4, France 



C. W. J. Beenakker 

Instituut-Lorentz, Leiden University, P.O. Box 9506, 2300 RA Leiden, The Netherlands 

(February 7, 2008) 



Abstract 



We calculate the joint probability distribution of the Wigner-Smith time- 
delay matrix Q = —ihS~^dS/de and the scattering matrix 5 for scattering 
from a chaotic cavity with ideal point contacts. Hereto we prove a conjecture 
by Wigner about the unitary invariance property of the distribution functional 
^[^(e)] of energy dependent scattering matrices S{e). The distribution of the 
inverse of the eigenvalues ti, . . . , tat of Q is found to be the Laguerre ensemble 
from random-matrix theory. The eigenvalue density p(r) is computed using 
the method of orthogonal polynomials. This general theory has applications 
to the thermopower, magnetoconductance, and capacitance of a quantum dot. 

PACS numbers: 05.45.-hb, 03.65.Nk, 42.25.Bs, 73.23.-b 
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I. INTRODUCTION 



The time delay in a quantum mechanical scattering problem is related to the energy 
derivative of the scattering matrix S{e). Although the notion of such a relationship goes back 
to the historical works of Eisenbud |l| and Wigner [0], the first to propose a matrix equation 
linking time delay and scattering matrix was Smith [Q . He introduced the hermitian matrix 
Q = —ihS~^dS/de, now known as the Wigner-Smith time-delay matrix, and interpreted its 
diagonal elements as the time delay of a wave packet incident upon the scatterer in one of 
N scattering channels. The eigenvalues ri, . . . , tat of Q are referred to as the "proper delay 
times" . 

In recent years, the issue of time delay and its connection to the scattering matrix has 
received considerable attention in the context of chaotic scattering [§-1^. Examples are 



the scattering of electromagnetic or sound waves from chaotic "billiards" or the transport of 
electrons through so-called "chaotic quantum dots" |2^. In both cases the scattering matrix 



S and the time-delay matrix Q have a sensitive and complicated dependence on microscopic 
parameters, such as the energy e, the shape of the cavity, or (in the electronic case) the 
magnetic field. As a result, a statistical approach is justified, in which one addresses the 
statistical distribution of S and Q for an ensemble of chaotic cavities, rather than the precise 
value of these matrices for a specific cavity. 

For the scattering matrix S, there is considerable theoretical and experimental 



25|^8| evidence that its distribution is uniform in the group of N x N unitary matrices. 



restricted only by fundamental symmetries, if the cavity is coupled to the outside world via 
ballistic point contacts. This universal distribution is Dyson's circular ensemble of random- 
matrix theory P^JSPj] . The fundamental problem is to calculate the distribution of the 
time-delay matrix Q for this case of ideal coupling. Known results include moments of tr Q 
P JT^JTB|] and the entire distribution for the case = 1 | ]I3| , p3 | of a scalar scattering matrix. 



In a recent paper we succeeded in computing the entire distribution for any di- 



mensionality A^ of the scattering matrix. We found that the eigenvalues ri , r2 . . . , r^v of Q 
are statistically independent of the scattering matrix S, with a distribution that is most 
conveniently expressed in terms of the inverse delay times 7„ = l/r„, 

P(7i, . . . ,7^) oc n l7. - 7,rn7r'e-/^--^/^ (1.1) 

i<j k 

Here (3 = 1,2,4 depending on the presence or absence of time- reversal and spin-rotation 
symmetry, th = 27r^/A is the Heisenberg time, and A is the mean level spacing of the closed 
cavity. The function P is zero if any one of the 7n's is negative. This distribution is known 



in random- matrix theory as the (generalized) Laguerre ensemble ||3^. The eigenvectors of 
Q are not independent of S, except in the case /3 = 2 of broken time-reversal symmetry. For 
any (3 the correlations are transformed away if Q is replaced by the symmetrized matrix 

Q, = -zhS-'/'^S~'/\ (1.2) 
oe 

which has the same eigenvalues as Q. The matrix of eigenvectors U which diagonalizes 
Qs = ?7diag(ri, . . . ,tn)W is independent of S and the r„'s, and uniformly distributed in 
the orthogonal, unitary, or symplectic group (for /3 = 1, 2, or 4, respectively). 
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In this paper we present a detailed and self-contained derivation of these results, focusing 
on those parts of the derivation that have not or only briefly been discussed in Ref. 



17 



We feel that a detailed presentation of this derivation is important, because it relies on 

(An 

The conjecture concerns the invariance of the 



an old conjecture by Wigner |3T[ that had not been proven before in the literature, 
unpublished proof is given in Ref 



00 



distribution functional P[S'(£:)] of the ensemble of energy dependent scattering matrices 
under unitary transformations S{e) — > VS{e)V' (with V, V two energy- independent unitary 
matrices). Our proof of Wigner's conjecture is based on the Hamiltonian approach to chaotic 
scattering which connects the scattering matrix S{e) of the open cavity with leads to 
the Hamiltonian 7i of the closed cavity without leads. It describes the so-called "zero- 
dimensional" or "ergodic" limit in which the time needed for ergodic exploration of the 
(open) cavity is much smaller than the typical escape time. In Ref. one of us used the 
Hamiltonian approach to derive the unitary invariance of the distribution function P{S) of 
the scattering matrix at a fixed energy. This unitary invariance is known as Dyson's circular 
ensemble 



29 



Wigner's conjecture is the much more far-reaching statement of the unitary 
invariance of the distribution functional P[S{e)]. 

The outline of this paper is as follows. In Sec. II we formulate Wigner's conjecture in 
its most general form for arbitrary dimensionality of the scattering matrix. (Wigner only 
considered the scalar case N = 1.) The proof follows in Sec. HI. In Sec. IV we show how 
the distribution of the matrices S and Q follows from this, now proven, conjecture. Then, 
in Sec. V, the density of proper delay times is computed using the theory of orthogonal 
polynomials. 

Historically, the Wigner-Smith matrix Q oc dS/ de was introduced to study the time evo- 
lution of a wave packet. This application is limited to the average delay time 0], essentially 
because knowledge of the full time dependence requires also higher derivatives of S than 
the first. Still, there exist physical observables that are entirely determined by Q. Some 
examples that have been discussed in the recent literature are the capacitance or admit- 
tance of a quantum dot [Q, and the thermopower [^. In addition to the energy derivative 
dS/de, one may study the derivative of 5* with respect to an arbitrary external parameter 
X. Such a parameter can represent the shape of the cavity, the magnetic field, or a lo- 
cal impurity potential. Recent examples of parametric derivatives range from conductance 
derivatives and charge pumping |]3^ in quantum dots to spontaneous emission [^] and 
photo-dissociation in optical cavities. The distribution of dS/dX can be obtained by 
methods similar to those used for the derivation of dS/de. A brief discussion is given in Sec. 



VI. We conclude in Sec. VII. 



II. PROBABILITY DISTRIBUTION OF THE SCATTERING MATRIX AND 

WIGNER'S CONJECTURE 

In order to describe the energy dependence of the scattering matrix S{6), we make use 
of the so-called Hamiltonian approach to chaotic scattering |^,33,40|. In this approach, the 



N X N scattering matrix S{e) is expressed in terms of an M x M random hermitian matrix 
Ti and a M x N coupling matrix W, 

S{e; H) = l- 2mW\e -H + i-nWW^Y^W. (2.1) 
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The hermitian matrix Ti. models the Hamiltonian of the closed cavity, its dimension M 
being taken to infinity at the end of the calculation. The coupling matrix W contains 
matrix elements between the scattering states in the leads and the states localized in the 
cavity. We distinguish three symmetry classes, labeled by the parameter (3: (3 = 1 (2) 
in the presence (absence) of time-reversal symmetry and /3 = 4 in the presence of both 
time-reversal symmetry and spin-orbit scattering. The elements of 7i are real (complex, 
quaternion) numbers for /3 = 1 (2, 4). The symmetries oiTi imply the following symmetries 
for S: S is unitary symmetric {(3 = 1), unitary (/5 = 2), or unitary self-dual (/3 = 4). 

For the problem of chaotic scattering one considers a statistical ensemble of cavities, 
which may be obtained by e.g. varying the shape of the cavity. We ask for the average 
of some function f[S{ei), . . . , S'(£:Ar)] over the ensemble. In the Hamiltonian approach, this 
ensemble average is represented by an integration over the matrix Ti, which is taken to be 
a random hermitian matrix with probability distribution -P(?^), 

if [Sis,], . . . , sisN)]) = I dnp{H)f[s{eu n),..., sisN, n)]. (2.2) 



Usually Ti. is taken from the Gaussian ensembles from Random Matrix Theory |30 



P(7^) cx exp ( - 5: F(E„) ) , nE) = ^^E' (2.3) 



n=l 



where Ei, E2, . . . , Em are the eigenvalues of Ti. The precise choice of F{E) is not important 
in the hmit M ^ 00 . 



The ensemble average defined by Eq. ( p.2| ) can also be formulated in terms of a distri- 
bution function of the scattering matrix S: 

if [3(6,), ^(£„)]) = j dSr... dSn Pie,, ^1, ... , S^,) /{S,, ...,8^,). (2.4) 

Here P{ei,e2, 5*1, ... , Sn) is the joint probability distribution of the scattering matrix 

5* at the energies ei, 62, ■ ■ ■ , Sn- The measure dS is the unique measure that is invariant 
under transformations 

S VSV , (2.5) 

where V and V are arbitrary unitary matrices {V' = V"^ for /? = 1 and V = for /3 = 4 
where T denotes the transpose and R the dual of a matrix). 

Statistical averages of the form (^) are completely characterized by the eigenvalues 



of W^W. For the case of ballistic point contacts (or "ideal leads"), these eigenvalues are all 
equal to MA/tt^ PJ2^J33|]. After application of suitable basis transformations on 5* and Ti 



one has therefore = 5^„(MA)^/^/7r for = 1, . . . , M, n = 1, . . . , N. In this case, the 
distribution P{e, S) of S{e) at one single energy is found to be particularly simple B,p4|, 

P{e] S) = constant. (2-6) 

Equation ( |2.6| ) is the starting point of the "scattering matrix approach" to chaotic scattering 
H . The ensemble of scattering matrices that is defined by the probability distribution 



4 



and the invariant measure dS is Dyson's circular ensemble from Random Matrix Theory 
29| . Depending on the symmetry class, we distinguish the Circular Orthogonal, Unitary, 



and Symplectic Ensembles, (COE, CUE, CSE), for /3 = 1, 2, and 4 respectively. 



The distribution function (|2.6|) is invariant under the transformation (|2.5| ). Wigner ||3T 
conjectured that the unitary invariance extends to the distribution functional P[S{e)], i.e. 
to the whole energy dependent ^-matrix ensemble. In other words, Wigner's conjecture is 
that (for any n) P{ei, 5*1, . . . , Sn) is invariant under a simultaneous transformation 

S.^VS.V, j = l,...,n, (2.7) 

where V and V are arbitrary unitary matrices {V = V"^ , for (3 = 1,4). 



III. PROOF OF WIGNER'S CONJECTURE 

Our proof of Wigner's conjecture consists of two parts. We first present an alternative 
random matrix approach for the energy dependence of the scattering matrix [O], and then 



show that in this approach the invariance property (|2.7|) is manifest. 

In the alternative random matrix approach, the role of the M x M matrix Ti is taken 
over by an M X M unitary matrix U, 

Sis) = f/n - f/i2 + U22y' f/21. (3.1) 

Here the matrices Uu, U12, U21, and U22 denote four subblocks of U, of size N x N, N x 
{M - N), {M - N) X N, and {M - N) x {M - N), respectively. 



U = ■ " ' , _ „ . (3.2) 




The energy-dependence of S enters through the phase factor exp(— 27ri£/MA). The matrix 
U is distributed according to the appropriate circular ensemble: U is distributed according 
to the circular ensemble, COE (CUE, CSE) for /3 = 1 (2,4). 

We now show that Eq. (|3.1| ) with M ^ iV is equivalent to the Hamiltonian approach 
( p.l|) . Equivalence of the circular ensemble and the Hamiltonian approach at energy e = 
was proven in Ref. [El]. This allows us to write the M x M unitary matrix U in terms of 



an M X M hermitian matrix Ti and an M x M coupling matrix W , 

U = l + 2'KiW\n-nxWW^Y^W. (3.3) 
The matrix Ti is distributed according to the Lorentzian ensemble [^], which in the large-M 



limit is equivalent to the Gaussian ensemble ( p.3| ). [The Lorentzian ensemble has F{E) = 
\{[3M + 2-/3) ln(M2 + vr^ E^).] 

The coupling matrix ly is a square matrix with elements W^n = 7r~^5^„(M A)^/^. We 
separate W into two rectangular blocks Wi and W21 of size M x N and M x [M — N), 
respectively, 

W={WuW2), (3.4) 
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stub 
5V^0 



FIG. 1. Physical motivation of Eq. (p.lj), explained in the text. 



and substitute Eq. ( |3.3| ) into Eq. (|3.1|). The result is 

S{e) = 1 - 2-KiWl{e -n-6H + inW^Wly^Wi 
where the M x M matrix 57i is defined as 



e \ }N 

e- (MA/7r)tan(7r£/MA) ) }M - N ' 



(3.5) 



(3.6) 



Apart from the term 57i, Eq. (|3.5| ) is the same as Eq. (|2.1| ) in the Hamiltonian approach. 
The extra term STi. is irrelevant in the limit M ^ cxd at fixed A^, because the diagonal matrix 
STi. contains only a finite number of matrix elements that are of order e, the others being 
of order e'^/MA. (One easily verifies that such a perturbation does not change eigenvalues 
and eigenvectors of Ti in the limit M ^ oo.) This proves equivalence of the alternative 
random- matrix approach ( p.l|) and the Hamiltonian approach (|2.1j ). 

It remains to derive the unitary invariance of P[S'(e)] from Eq. ( |3.1| ). As the matrix U is 
chosen from the circular ensemble, its distribution is invariant under a mapping U — U'UU", 
where U' and U" are arbitrary M x M unitary matrices {U"U'* = 1 if /? = 1,4). We now 
choose 



U' 




U" 




(3.7) 



where V and V' are arbitrary unitary N x N matrices {y*V' = 1 if /3 = 1,4). Invariance 
of P{U) under the transformation U U'UU" with the choice (|3.7| ) implies, in view of Eq. 
( ^.1| ), invariance of P[S'(e)] under the transformation 



S{e) VS{e)V', 



(3.8) 



as advertized. 

A central role in our proof of Wigner's conjecture is played by the random-matrix model 
( p. ID , which expresses the energy-dependent N x N scattering matrix S{e) in terms of an 
energy-independent M x M random unitary matrix U . This model was first proposed in 
Ref. [|l6l, following a more physical reasoning than the formal derivation given above. It is 
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insightful to briefly repeat the reasoning of Ref. |T6|. We consider the scattering matrix S 
at energy e = and then shift the energy by an amount e. Equivalently, we can replace 
the energy shift e by a uniform decrease 6V = —e/e of the potential V in the quantum dot. 
Since the quantum dot is chaotic, we may localize 6V in a closed lead (a stub), see Fig. |l|. 
The stub contains M — N ^ N modes to ensure that it faithfully represents a spatially 
homogeneous potential drop SV. The scattering properties of the system consisting of the 
dot, the A^-mode lead, and the (M— A^)-mode stub are described by the (M— A^)-dimensional 
reflection matrix of the stub and the M-dimensional scattering matrix U of the cavity 
at energy e = 0, with the stub replaced by an open lead. The advantage of localizing the 
potential shift 6V inside the stub is that the energy-dependence of is very simple, 

r,{e) = _e-2--/A^^. (3.9) 

The matrix U is energy independent by construction. It is taken from the circular ensemble, 
because scattering from the dot is chaotic. (The matrix is not random.) Expressing the 
scattering matrix S{e) in terms of U and rs{e) then yields Eq. (p.l|) . 



IV. DISTRIBUTION OF THE TIME-DELAY MATRIX 

In this section we combine Wigner's conjecture for the invariance properties of the scat- 
tering matrix ensemble and the Hamiltonian approach to chaotic scattering, to compute 
the joint distribution P{S, Q^) of the scattering matrix S{e) and its symmetrized energy 
derivative 

= -,ns~"^^s-''\ (4.1) 

OS 

The matrix is a real symmetric (hermitian, quaternion self-dual) matrix for = 1 (2, 
4). The joint distribution P{S, Q^) of S and is defined through a relation similar to Eqs. 



J) and dJ) 



{f[S{e),Q,{e)]) = J dnP{n)f[S{e,n),Qe{e,H)] 

^ J dSdQJ{S,Q,)P{S,Q,), (4.2) 

where /(S*, Qe) is an arbitrary function of S and Qg. 

Application of Wigner's conjecture is the key to the calculation of P{S, Qg). To see this, 
we first consider the distribution function P{e2 — £i; 5*1, 5*2) of the scattering matrix S{e) at 
the energies ei and £2, for which the invariance property ( |2.7|) implies 



P(£2 - ^i; 5i, ^2) = P(£2 - £1; VS^^', VS2V') 

= P{e, - 6,; -1, -Sr'/' S2 Sr'/'). (4.3) 

— 1/2 

For the second equality we have chosen V = V = i Si . In the limit £2 ~^ ^1 = ^, one 
finds from Eq. ( [4.3|) that the joint distribution P{S, Qg) of the scattering matrix S{e) and 
the symmetrized time-delay matrix Qg does not depend on the scattering matrix S, 
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P{S,Qs)=Pi-l,Qs). (4.4) 

Hence, to find P{S, Q^) it is sufficient to calculate the integral ( |4.2| ) for a function /(S*, Q^) 
of the form 

/(5,Q,) = 5(5, -l)/(g,), (4.5) 

where f{Qe) is an arbitrary function of Qg and 5(5*, —1) is the delta function at S' = — 1 on 
the manifold of unitary symmetric (unitary, unitary self-dual) matrices for (3 = 1 (2, 4). 
We now turn to the evaluation of the integral ( |4.2|) , with / of the form (|4.5| ), using the 



Hamiltonian approach. We note that (|2.1| ) can be rewritten as 



5 = -l + -^— , K = 7:W^^-W. (4.6) 
1 + iK e —ri 



The energy derivative of S is given by 



^ = 2m i W^- — ^-—^W ^ . (4.7) 



We decompose the hermitian matrix TC in Eq. ( [4.2|) into its eigenvalues and eigenvectors: 
Ti. = ipEtp^ where E is the diagonal matrix containing the eigenvalues Ep of Ti as entries 
and ip is the orthogonal (unitary, symplectic) matrix of the eigenvectors of 7^ for /3 = 1 (2, 

4), 

(/(5,Q,)) = J dtp J dE, . . .dEM Jm{E)P{E) f[S{e;tljEtlj^),Qg{e;tljE^^)], 

M 

Jm{E)= ni^p-^.l^- (4.8) 

P<(T 

Here dip is the invariant measure on the orthogonal (unitary, symplectic) group and Jm{E) 
is the jacobian for the variable transformation Ti ^ ip^E [Q. 

We now make the special choice ( |4.5| ) for f{S,Qs). In terms of the matrices ip and E, 
the matrix K from Eq. (|4.6| ) reads 

;.„„.M^^^. (4.9) 

^ p=i ^-Ep 

Inspecting ( [4.6| ), we see that the limit 5^—1 corresponds to the case when all N eigenvalues 
of K diverge or, equivalently, when at least out of M eigenvalues Ep of Ti tend to the 
energy e. Let us label these eigenvalues Ei, . . . , E^, i.e. \Ep — e\ — in the limit S —1 
for p = 1, . . . , N. In this limit, the sum in (|4.9| ) can be restricted to the first N terms and 
both S and dS/de do not depend on the eigenvalues Ep and on the matrix elements tpmp 
with m > N OT p > N. Therefore we can perform the integration with respect to the latter, 
resulting in 

{S{S,-l)f{Qg)) = J dQJ{Qg)P{-l,Qg) 

= Jd^P{^) JdE,...dEr,P{E,,...,Er,)S{S,~l)f{Qg), (4.10) 
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where \1/ is the N x N matrix containing the rescaled eigenvector elements tpmpM^I'^ for 
1 < m, p < N , the integration measure d^! = nm,p=i ^^mp' -f(^) is the distribution of 
\t' after integration over the remaining matrix elements of ip, and P{Ei, . . . , E^) is the 
distribution of the eigenvalues Ei,. . . ,En after integration over the remaining eigenvalues 
Ep, p> N. Near S = -1 we have from Eqs. (|;6D, Q, and (U), 

^ = -1 + (27rz/A)^^""Miag {Ei - e, . . . , En - (4.11) 
Q, = r^vl/t-ivp-i, (4.12) 

where th = 2nh/A is the Heisenberg time. In the limit M ^ N, the matrix elements of ^ 



are Gaussian distributed with zero mean and unit variance 30 



P(^) oc exp[-(/3/2)tr ^^^]. (4.13) 

Because of the delta function 6{S, —1), it is sufficient to know P{Ei, . . . , Ej^) for \Ep—e\ 0, 
p = 1,...,N. In this case, P{Ei,...,En) is entirely determined by the Jacobian 
Jn{Ei, . . . ,En), up to a numerical factor which results from the integration over the re- 
maining M — N eigenvalues of Ti, 

N 

P{E,,...,EN)^Y[\Ep-E'^f ■ (4-14) 

p<cr 

We have now succeeded in transforming the integral ( |4.2| ) over an M x M matrix Ti to 
the integral ( [4.1CI| ), which is formulated in terms of N x N matrices only. At this point, 



it is convenient to rewrite the integration over the eigenvalues Ei, . . . ,En as a matrix 
integration to avoid the singularity of P{Ei, . . . , En) at Ei = . . . = En = £■ Hereto we 
exploit that the distribution of \1/ is invariant with respect to a transformation of the type 
\I' —>■ O"^, O being an orthogonal (unitary, symplectic) matrix for (3 = 1 (2, 4). We apply 
this transformation for arbitrary O and average over O with respect to the invariant measure 
dO. Then, we substitute the hermitian N x A^-matrix H = OMiag {Ei, . . . , EN)0—e and the 
corresponding Jacobian exactly cancels the singular level-repulsion factor of P{Ei, . . . , En), 

(5(5,-l)/(Q,)) oc I d^Pi^) I dH5{S,-l)f{Qe) 

oc j d^P{^) J dH6{^^~^ H^-^)f[rH'^^'^'^'^]. (4.15) 

(The proportionality sign indicates that we have omitted a normalization constant.) Three 
more variable transformations are needed to calculate the integral (|4.15|) . First, we replace 
the integration variable H hj H' = if The Jacobian det(\I^ \I'^)(^'^+^"'^)/^ of this 
transformation can be derived using the singular value decomposition of \E' |^2|]. The integral 
( [4.15| ) becomes 

{6{S,-l)f{Q,)) oc Jd^P{^) JdH'6{H') det(^^t)(^^+2-W2j[^^$t-i^-i] 

= J rf^e-(^/2)tr*vI-t^g^(^^t)(iV/3+2-;3)/2^[^^^t-l^-l]^ ^4_^Q) 
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Next, we note that the integrand depends on the hermitian matrix T = "^"^yrn only, and 
replace the integration variable \1/ by F = \l/\l/yr/^. The Jacobian of this transformation 
provides an extra factor oc det(r)^/^~^^^(r), where 9{T) = 1 if all eigenvalues of F 
are positive and 0{r) = elsewise. Finally, we replace F by its inverse Q = F~^. Since this 
variable transformation has a Jacobian det((5)^^"*'^~^, we arrive at 

{6{S, -l)fm) oc / dQe{Q) det(g)-3^^/2-2+^e-^^-*=-«-'/V(Q). (4-17) 

Using Eqs. ( |4.4| ) and ( |4.1CI|) we thus find the joint probability distribution P{S, Qe) of the 
scattering matrix S and the symmetrized time delay matrix Q^, 

P(S,g,) =P(-1,Q,) =^^(g,)det(Q,)-3^^/2"2+/3e-'3-«trQrV2_ (418) 

The corresponding distribution of the eigenvalues ti, ... ,tn of Qs, the proper delay times 
then reads 

n<m n 

Alternatively, the distribution of the rates 7^ = l/r,- is given by 

P(7l,...,7^) « n \^n-lraf ViQ{ln)lt'"^ e-'^--W2 . (4.2O) 
n<m n 

This distribution is known in random-matrix theory as the generalized Laguerre ensemble 
1 3^] . The eigenvectors of Qs are uniformly distributed according to the invariant measure of 
the orthogonal (unitary, symplectic) group. 



V. DENSITY OF PROPER DELAY TIMES 



The distributions ( [4. 191 ) and ( [4. 20] ) have a form which allows to apply the standard 



technique of orthogonal polynomials to determine the density or the correlation functions 
of the eigenvalues of (or F = Q^^). The A^-dependent exponent is somewhat unusual 
and implies that the set of orthogonal polynomials is different for each value of A^. Such 
ensembles have been studied in mathematical physics [^5i^ |. 



We restrict ourselves to the simplest case /5 = 2 of broken time reversal symmetry. For 
a simplified notation, we use the dimensionless escape rates Xn = 'JuTh = 'Tn/'Tn which are 
distributed according to 

p{xi,. . . ,xn) \{{xn- x^f \{wN{Xn) , wn{x) = x^ e'"" . (5.1) 

n<m n 

The generalized Laguerre polynomials (x) are orthogonal with respect to the weight 
function wn{x). The method of orthogonal polynomials relates the correlation functions of 
the x„'s to these polynomials. Here we only consider the density 

p{x) = I dxi . . . dxN pixi, . . . , Xn) 5{x — Xn), (5.2) 

•^0 n=l 
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0.5 1 1.5 2 
Nt/t„ 

FIG. 2. Density of proper delay times. Shown are the densities for N = 1 (full line), = 3 
(dashed line), both computed from Eq. ([5.3|), and the limit N ^ oo (dotted line) given by Eq. 



which is given by the series 



n=0 



(N + n)] 



Using the recurrence relation 



(n + 1) Li%\{x) = {2n + N + l-x) L^^\x) -{n + N) L^^\{x) 



{N) 



(5.3) 



(5.4) 



with Lq (x) = 1 and l[^\x) = + 1 — x, it is possible to evaluate Eq. ( |0| ) efficiently for 
small A^. For large N, one can use an asymptotic expansion of the Laguerre polynomials to 
find the closed expression 



p{x) 



271X 



V-N^ + QNx - x^ . 



The corresponding density for the proper delay times t = th/x thus reads 

N n TH, 



Pir) 



27rr2 



r^-r)ir-r^) , r± = ^(3±V8). 



(5.5) 



(5.6) 



The limit — oo is rapidly approached for A^ ^ 3, see Fig. 2. 

It is instructive to compare the density of the proper delay times r„ with the density 
of partial delay times studied by Fyodorov, Sommers, and Savin [|,|l3. The partial delay 
times tn are defined in terms of the eigenvalues e*"^" of the scattering matrix S, 
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0.5 1 1.5 2 

Nt/T„ 

FIG. 3. Density of partial delay times, from Ref. |l^. Shown are the densities for = 1 (full 
line), = 3 (dashed line) and N = 100 (dotted line). The limit N ^ 1 corresponds to a gaussian 
with mean 1 and width l/^ 




h 



d(t)n 

de 



(5.7) 



Compared to the definition of proper delay times as eigenvalues of the Wigner-Smith matrix, 
we see that for the partial delay times the order of the operations of energy derivative and 
diagonalization is reversed. This order is irrelevant for the sum J2n '^n = J2n tn, corresponding 
to the first moment of the density, but higher moments are different (unless N = 1). The 
density of partial delay times was obtained in Ref. @Jl^ by the super symmetric approach, 



Pit) 



{N-iy.TH 



t 



N+2 



-TH/t 



(5.8) 



Notice that this density corresponds to the contribution from the first Laguerre polynomial 
in the summation ( ^.3[ ). For N >2, there are also contributions from higher order Laguerre 
polynomials and therefore the two densities do not coincide. In Fig. |^, the density of the 
partial delay times is shown for = 1, 3, 100. For ^ 1, Eq. ( |5.8|) becomes a gaussian of 
mean 1 and a width l/\fN. The qualitative difference between p(r) and p{t) is due to the 
absence of level repulsion for the partial delay times. Unlike the proper delay times, we do 
not know of a physical quantity that is determined by the partial delay times. 



VI. GENERAL PARAMETRIC DERIVATIVES 

So far we have restricted our discussion to the energy derivative of the scattering matrix. 
It is also of interest to know the derivatives with respect to external parameters. Examples 
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of such parameters are a magnetic field, tlie sliape of the cavity, or a local potential. In 
this section we present and prove an extension of Wigner's conjecture to the ensemble of 
scattering matrices S{e, X) that depend both on energy and on an external parameter X, 
and then use it to compute the distribution of the symmetrized derivative 

Qx = -iS~'/'^S-"\ (6.1) 

In the Hamiltonian approach, the parameter dependence of the scattering matrix is 
modeled through Eq. ( p.l|) with a parameter dependent random hermitian matrix ?-^(X), 

S{e, X-n) = l- 27TiW^[e - n{X) + mWW^Y^W, 

n{x) = n + xn'. (6.2) 

The matrix 7i' is a random hermitian, Gaussian distributed, symmetric (antisymmetric) 
matrix if X denotes a shape change (magnetic field), or Ti'^j = 6ir6jr if X denotes the 
local potential at site r. (In random-matrix theory, the site r is represented by an index 
1 < r < M, see e.g. Ref. ^8[|.) As we argue below, the invariance property ( |2.7|) , which 
was conjectured by Wigner for the ensemble of energy-dependent scattering matrices, is also 
valid for the more general ensemble of parameter- and energy-dependent scattering matrices 
S{e,X). That is to say, the distribution functional P[S{6, X)] of the ensemble of matrix 
valued functions S{6,X) is invariant under transformations 

S{e,X) ^VS{e,X)V', (6.3) 

where V and V are arbitrary unitary matrices {V*V' = 1 if [3 = 1,4). The proof of 
this invariance property goes along similar lines as in Sec. |T|: First, one shows that the 
Hamiltonian approach (|6.2|) for the ensemble of parameter dependent scattering matrices is 



equivalent to a formulation |48| in terms of unitary matrices of a form similar to Eq. ( |3.1|) 



S{e,X) = t/n + t/i2 [e-2-(-^^')A^^ - U,,. (6.4) 



Second, one verifies that in the formulation ( |6.4|) the invariance property (|6.3|) is manifest, 
completing the proof. 

Using the invariance property (|6.3|), we can now compute the distribution of Qx- For 
simplicity, we restrict ourselves to the case that the matrices Ti and Ti' have the same, 
Gaussian, distribution. As in the case of the energy derivative of the scattering matrix, see 



Sec. IV, it is sufficient to consider the special point S = —1. From Eqs. (WM) and (|6.2|) we 



find that at 5* = — 1 



= ^fJ^-'H'^-\ P{H') oc e-^*^^"/^^ (6.5a) 

H'^n = -{rHim-"^Y.rm,Ku^n.- (6.5b) 



Repeating the steps outlined in the previous section, we find that the distribution of Qx is 
Gaussian with a width set by Q^, 



13 



P{S,Qe,Qx) oc (detQe) 



-2/3Af-3+3/3/2 



exp 



tr 

2 



-1 

e 



(6.6) 



The reason why the time-delay matrix sets the scale for the matrix Qx characterizing the 
response to the external paramater X can be understood in a picture of classical trajectories 
49| : For long delay times, the scattering properties are more senstive to a perturbation of 



the system than for short delay times. Such an explanation in terms of classical trajectories 
is valid in the semiclassical limit iV — oo only. Our exact result (3^) makes the relation 
between time delay and parameter response precise in the fully quantum mechanical regime 
of a small number of channels N . 

The invariance property (|6.3|) and the Gaussian distribution ( |6.6|) also hold if the scatter- 
ing matrix depends on more than one external parameter. If the scattering matrix depends 
on parameters Xi,. .., the matrices Qx^ (j = l,...,n) all have the same Gaussian 
distribution, with a width set by the symmetrized time-delay matrix Q^, 



P{S, Qe, Qxi, • • • , QxJ oc (det Q, 



X exp 



■,-N/2-il3N+2-/3)in+2)/2 



B 1 " 

-ftr \rHQ:' + -Y:yHQ:'Qx, 



(6.7) 



VII. CONCLUSION 

In summary, we have presented a detailed derivation of the joint distribution P(5', Q^) of 
the scattering matrix S and the symmetrized time delay matrix = —ihS~^^^{dS/de)S~^^'^ 
of a chaotic cavity with ideal leads. Our result, which was first reported in Ref. [|1^], reads 

P{S, Qe) OC eiQe) det(ge)-=^'3^/2-2+/3 ^-f^r^tr Q-^ /2 ^ (7^^) 

where th = 2iTh/A is the Heisenberg time and A is the mean level spacing in the cavity. The 
distribution P{S, Qe) depends on the eigenvalues Ti, . . . ,T]y of Qe only; it does not depend on 
S, nor on the eigenvectors of Qe- Our derivation was based on an old conjecture by Wigner 
31| , which we have proven in this paper, that the distribution functional P[S'(£)] of the 



ensemble of energy-dependent scattering matrices is invariant under unitary transformations. 
We generalized Wigner's conjecture to the dependence of S on an external parameter X and 
derived the distribution of parametric derivatives. 

Throughout this paper we assumed ideal coupling of the chaotic cavity to the electron 
reservoirs. This requires ballistic point contacts. For the case of non-ideal coupling, i.e. if 
the point contacts contain a tunnel barrier, the distribution of the scattering matrix S is 
given by the more general Poisson kernel instead of the circular ensemble For this 



situation, application of Wigner's conjecture is not correct and our method can not be used 
to determine P{S,Qe)- A possibility to overcome this difficulty is to separate the direct 
and the fully ergodic scattering processes, by expressing the scattering matrix for non-ideal 
coupling as a composition of an x scattering matrix Serg of a fully ergodic and chaotic 
scatterer and a scattering matrix S'dir of dimension 2N x 2N describing the direct scattering 
processes [P3|,H2|,ET| , 
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S = r + t\l - S,rgr')-^Sergt, (7.2) 
where the matrices r, r', t, and t' are N x N submatrices of S'dir, 



-^dir = • (7.3) 




The energy dependence of S'dir can be neglected on energy scales corresponding to the inverse 
dwell time inside the cavity. Thus one obtains a relation between the statistical distributions 
of the time delay matrices for the cases of non-ideal and ideal coupling. What remains is the 
cumbersome mathematical problem to perform the corresponding transformation of the dif- 
ferent matrix variables. To our knowledge, no closed-form formula for the joint distribution 



of S and Qe is known, except for the case = 1 ||5T[ . 

The supersymmetric approach used in Refs. P JT^JT^J53| ] is typically more general and 
systematic concerning the issue of non-ideal coupling, since there the particular symmetries 
of the scattering matrix ensemble in the ideal coupling case are not directly exploited. On the 
other hand, this technique is usually restricted to the computation of certain Green function 
averages only, and does not give access to the full distibution of the time-delay matrix. We 
believe that, in principle, in the supersymmetric approach it should be possible to determine 
averages of the form (tr (Q")) for arbitrary n (and hence the density of delay times) and for 
arbitrary coupling, without having to increase the dimension of the supermatrices involved. 
However, in order to get access to the full distribution of the matrix Q^, or at least to 
moments of the type ((trQ)"), one has to increase the dimension of the supermatrices 
accordingly, which makes it apparently impossible to go beyond the case n = 2 0,|12|,|13|. It 



is amusing that the (now proven) half-a-century old conjecture of Wigner provides a route 
to the simple result ( |7.1| ) that seems unreachable by later supersymmetric techniques. 
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